Copyright © 1972 American Telephone and Telegraph Company 

The Bell System Technical Journal 

Vol. SI, No. 6, July-August, 1972 

Printed in U.S.A. 



Queues With Mixed Renewal 
and Poisson Inputs 

By A. KUCZURA 
(Manuscript received February 2, 1972) 

In a queueing system with two independent in-put streams, as exists, 
for example, when first-routed and overflow traffic streams are offered 
to a common sender-group, the state of the system encountered by the two 
different types of customers upon their arrival will generally be different. 
Consequently, in a system where delayed customers wait for service, the 
service rendered to the individual streams may also be different. 

The delay distribution in a single-server queue for each type of customer 
is derived under the assumption that one stream is Poissonian and the 
other is described by a renewal process. The difference in service received 
by the two streams is examined with the aid of numerical examples for 
two interarrival time distributions of the renewal stream. We show for 
two cases that a practical indicator of service received by the renewal cus- 
tomers is the coefficient of variation of their interarrival time distribution. 
If the coefficient is less than unity, then the renewal customers receive 
better service than the Poissonian customers. The converse is true when the 
coefficient exceeds unity. 

The stationary distribution of the number of busy servers in an infinite- 
server system as seen by the two types of customers is also derived. 

I. INTRODUCTION 

The concept of a piecewise Markov process 1 is used to analyze two 
queueing systems the inputs of which are composed of two independent 
streams. One of the streams is Poissonian with intensity X and the other 
(called a GI stream because of its General Independent Distribution 
of intervals between arrivals). is assumed to be a renewal process with 
intensity v. We assume the service times of all the customers are in- 
dependent and identically distributed according to an exponential 
distribution with mean n~ l . Such models are denoted by GI + M/M/c 
in Kendall's notation, where the "M" refers to the Markovian character 
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of Poissonian arrivals and of exponential service and c refers to the 
number of servers. 

The state of the system (number of customers waiting and in service) 
seen by an arriving Poissonian customer will generally differ from that 
of the GI customer. Consequently, in a system where delayed customers 
wait for service such as in the GI + M/M/l queue, the service received 
by the two types of customers will differ. Whether the GI customers 
receive better or worse service than the Poissonian customers depends 
on the variability of the interarrival times of the GI stream. 

In Section II, we analyze the GI -f M/M/l queue. The delay distri- 
bution with order-of-arrival service for the two types of customers is 
derived for GI streams the interarrival time distributions of which 
have rational Laplace-Stieltjes transforms. The coefficient of variation 
of the interarrival time distribution is introduced as a practical measure 
of the variability properties of a stream of customers. Its usefulness 
in predicting service is evaluated through some numerical examples 
and its relation to the common measure called peakedness is examined. 
In Section III, in order to describe the intrinsic character of the two 
streams of customers as they would be observed in a system without 
delay, we derive the stationary distribution of the number of busy 
servers seen by the two types of arriving customers in a GI + M/M/ oo 
system. 

II. THE GI + M/M/l QUEUE 

Let Y(t) be the number of customers in the system (those waiting 
and in service) at time t. Since the GI stream is a renewal process and 
Y(t) is Markovian between any two consecutive arrival epochs of the 
GI customers, \Y(t), t ^ 0} is a piecewise Markov process 1 with state 
space {0, 1, 2, • • • }. The regeneration points are the arrival epochs 
of the GI customers. Thus the distribution of the length of the 
Markovian segments is given by A(£), the interarrival time distribution 
of the GI customers. The regeneration matrix is given by 



<P«) = 



10 
10 
1 



(1) 



The elements of this matrix are the transition probabilities across a 
regeneration point; that is, p,, is the probability that immediately 
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after regeneration the process is in state j, given that, immediately 
prior to regeneration, the process was in state i. 

The Markov process operating within the segments is a birth-death 
process with birth rate X and death rate p, the same for all segments. 
The transition probability functions of this process are given in Ref. 
2, page 13. 

PM - P i "'"'V ,u " ,, {/ 1 , 1 (2<\/v) + P -*i, +l+1 (2*Vv) 

+ (1 - p) E p* ( * +1 7, + , +4+1 (2« Vxp")} , 

where p = X/p (X/p < 1) and /,(£) is the modified Bessel function 

The regeneration matrix (1), the distribution A(£), and the transition 
functions P,,(0 determine the piecewise Markov process completely. 
We will first derive the distribution of the state of the system seen by 
an arbitrary GI arrival and then use the rate conservation principle 
to find the stationary distribution 

q { = limP(F(0 = j\ 7(0) = »j, 

j = 0,1,2, ••• , for all i, (2) 

which is the same as the distribution of the state of the system seen 
by a Poissonian arrival. Having found the two distributions, we can 
readily determine the individual delay distributions for an order-of- 
arrival service discipline. 

2.1 Delay Sustained by GI Customers 

Let ( Pi ■} be the stationary distribution of the Markov chain imbedded 
at points immediately preceding a GI arrival. If r,, is the one-step 
transition probability from state i to state j, then 

Tu = r Pi+1,,® dA(Q, i,j = 0, 1,2, ••• . 
•'o 

The distribution {p,j satisfies the Chapman-Kolmogorov equations 
Pi = Y,Virn , i = 0, 1, 2, ••• , 



1308 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1972 

and hence the solution to the linear system 

Vi = f S P«P«+i.i® dA(Q, j - 0, 1, 2, • • • , (3) 

Jo t=0 

subject to the normalization condition 

CO 

1=0 

is the stationary distribution [p,-\. 
Let 

oo 

p(«) = EpX 

be the probability generating function of jp,}. Multiplying (3) by u' 
and summing over all j, we obtain 

p(«)= P I>.r, +1 ( M , *) eM®, (4) 

^O t=0 

where 

r.(«,*) = EPottK 

1=0 

is the probability generating function of P.- ,-(£), j = 0, 1, 2, • • • . This 
function is not readily available. Its Laplace transform 



7 ,-(w, s) = f e-'r.Cu, <£, 
Jo 



though, has the following form 2 : 

u i+1 - (1 - u)y i+ \s)/[l - ,(«)] 



where 



and 



7,(m, s) = . , , Vl 

u[s — h(u)\ 



h(u) = - (1 — u)(fi — Xw), 
u 



,s X + m + g - V(X + m + s) 2 - 4X M 

We will transform the real integral 

f r,- +1 ( W , Q dA{$) (5) 

Jo 
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appearing in (4) into a complex integral involving 7, + 1 (w, s) and a(s), 
the Laplace-Stieltjes transform of A(£). 
If we set 

dA(0 = e~ ai dB(0, a > 0, 

then (5) becomes the Laplace-Stieltjes transform integral 

f e- a tr i+l (u,Z)dB®. (6) 

Jo 

But the transform of a product is the complex convolution of the 
transforms of its factors and so (6) becomes 

1 r c+,w 
2^ J y i+ i(u, z)0(a — z) dz, 

where c is any positive number, i = y/ —\, and 

ft*) = fe-tdB®, (7) 

the Laplace-Stieltjes transform of B(£). But 

/3(s) - f° e-^- a) dA(0 = a(s-a), 
Jo 

and we finally have 

£° T i+1 (u, 5) dA(Q = ~ £2 Tf+i(«. «>><-«) &• (8) 

Substituting for y i + x {u, z) on the right-hand side of (8), and using the 
resulting identity in (4), we obtain the following integral equation 
for p(u): 

Piu) ~ 7T~- / r — ,/ \i \<*(—z) dz, (9) 

2-rn J c - iK L u[z — h(u)] J 

where 

V\Z)P[V(Z)] 



H(z) = 



1 - V (z) 



We will evaluate the complex integral in (9) for the class of a(s), the 
members of which are rational functions; but first we note some prop- 
erties of the integrand which suggest a contour to be used in applying 
the calculus of residues. 

Since P ,-,-(£) is a probability function, T,(w, £) is uniformly convergent 
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for | u | ^ 1 and 7,(w, z) is holomorphic for | u | ^ 1 and ®.(z) > 0, 
where (R(z) denotes the real part of z. Hence the bracketed part of the 
integrand in (9) is holomorphic f or | u | ^ 1 and ®.(z) > since it is 
the convergent series 

CO 

Sp,-7. + i(w, z). 

Since A{£) is a probability distribution function, a(z) is holomorphic 
for 01(2) > 0. For (R(z) < 0, a (z) may or may not be holomorphic. 
The predominant case is where a(z) is meromorphic in the half-plane 
(R(z) < and we shall address ourselves to this case. 

Let — Zi , — z 2 , • • • , — z n be the poles of a(z). Since the poles of a(z) 
are in the left-half plane, the poles of a (—2) are in the right-half plane. 
Hence, the integrand in (9) is meromorphic in the right-half plane 
and we can use the residue theorem to evaluate the integrand over 
the contour consisting of the line (c + iR, c — iR) and a semicircle of 
radius R in the right-half plane which connects c — iR with c + iR. 
We choose c and R such that all the poles of a ( — 2) are interior to this 
contour. 

Since a(z) is meromorphic, we can write 



, v Q.C0 



*n ' 



(2 + 2l )*'(2 + Z*)"* •■■ (Z+Zn) 

where fc, is the order of the pole at 2, , Q m (z) is a polynomial of degree m, 
and we assume that m + 1 ^ fc, + k 2 + • • ■ + k n . Two examples are: 
(i) A (£) is the gamma distribution with density 

«©- (jd^ (fe * )fc " V *'*' * >0 ' 
and transform 

■« = fe)' ; . C10) 

(m) ^4 (^) is the mixture of exponentials (hyperexponential) 
A(£) = E «.(1 - e""*), * > 0, 

i = l 

where a, > 0, r, > and o 1 + a a +---+a B = l, and transform 

^-S^ft- (11) 
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If we now let R — > oo , then the contour integral along the line z = c 
tends to the desired limits and it can be shown that the integral along 
the semicircle tends to zero. Hence, using Cauchy's integral formula 
with 

4i m \ - tt 2 pfc) ~ (1 ~ u)H{z) n . . 
m " u[z- h(u)] W " C Z) ' 



we have 

J_ p' w (-1)%) 

2ri h-i« (z - Zl ) k '(z - z 2 f> ■ ■ ■ (z - z n ) 



dz 



^-ir'E^T^-fe), (12) 



where k = fc, + k 2 + ■ ■ • + K and 

m 



gM = 



{z _ Zi y> . . . {z _ z^y-^z - ,. +1 )*' + ' ■ • • (z - z n ) 



If all the poles are simple, such as in example {it), then the integral 
in (12) is equal to 

(-D n+1 EfcOO. (i3) 

J = l 

If z x is the only pole and it is of order k x , such as in example (i), then 
the integral in (12) is equal to 



(*i " I)' 



f"-"W. (14) 



Detailed analysis for the case k = 2 will be carried out later when we 
give numerical examples. The case of two simple poles will also be 
analyzed. 

In general, after the integral has been evaluated, and the result 
substituted into (9), an equation in p{u) results. This equation can be 
solved for p(u). The values of the unknown function H(z) at the poles 
2, and the values of its derivatives may be determined by first applying 
the normalization condition ;)(1) = 1 and then, since p(u) is holo- 
morphic in | u \ ^ 1, forcing the zeros of the numerator in the unit 
circle to coincide with the zeros of the denominator in the expression 
for p{u). This procedure will be illustrated later by an example. 

We return to the delay distributions. Let W be the delay, or waiting 
time, from the arrival instant until the beginning of service for an 
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arbitrary GI customer. The GI delay distribution is defined by 

D(t) = P[W ^ t], t^ 0. 

If D t (t) is the conditional delay distribution, given that the arriving 
customer finds ;' other customers in the system, then 

D(t) = £ Z>,(i)p, • 

1-0 

Obviously, D (t) = 1. Since service is in order of arrival, 

and Dj(t) for j > 1 is the convolution of j identical exponential distri- 
butions, each with mean jT\ If 5(s) is the Laplace-Stieltjes transform 
of D(t), then 



i-l Jo 



i -I 

since the transform of a convolution is the product of the transforms. 
Consequently we have 



•» - fe) • 



This equation can now be inverted to obtain the delay distribution 
for the GI customers. 

2.2 Delay Sustained by Poissonian Customers 

Let E(t) be the delay distribution of Poissonian customers; that is, 
E(t) is the probability that an arbitrary, arriving Poissonian customer 
will be delayed no more than t units of time. If e(s) is its Laplace- 
Stieltjes transform, then, by the same argument as above, 

where 

00 

g(«) = Y, Qiu' 

1=0 

is the probability generating function of |g,}, the distribution of state 
seen by an arbitrary, arriving Poissonian customer. 
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Using the rate conservation principle 1 we can write 

vpj + Xg y = nq i+1 , j = 0, 1, 2, • • • , (15) 

where the left-hand side is the asymptotic rate of transition from j 
to j + 1 and the right-hand side is the asymptotic rate of transition 
from ;' + 1 to j. Multiplying both sides of (15) by u' and summing over 
all j, we obtain the following relation between p(u) and q(u): 

vp(u) + \q(u) = ^ lq(u) - go]. 
Thus, if p(u) has been found, then q(u) is given by 

q{u) = ft + t-rtu) _ 

where a = vj y. and p = X/p. Applying the normalization condition 
g(l) = 1, we get q , i.e., 

q n = 1 — p — cr, 

and hence 

«M = i - e =-; + V£M . (16) 

J. pw 

In terms of the probabilities themselves, we solve (15) for \q t \ and 
obtain 

q, = (1 - p - a)p j + a 2 P.P'" ,_I , | - 1, 2, ... . 

i-0 

2.3 Some Relations Between the Two Results 

If Z, is the mean of the distribution [p f ] and Z 2 is the mean of the 
distribution [q,], then using (16) and differentiating, we obtain 

fc - r^7 + r^7 (1 + w- (17) 

Note that the first term on the right-hand side of (17) is the mean 
number of customers in the system in the M/M/l queue with traffic 
intensity p. Hence, the second term may be considered to be an increase 
in the mean number due to the presence of the GI customers. 

Since the probability distribution functions p(u) and q(u) are related 
through (16), the transforms of the delay distributions are related 
through the following equation: 
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«0 - 7T7^ (1 - " - a) + *~+7^ 8(8) - (18) 

Taking the inverse transform of both sides of this equation and re- 
arranging the terms, we obtain 



(19) 



E{t) = 1 - p.-*-* 1 - all + m £ [P - J>(« " 0]e- (,, - w * #}• 

If n is the mean delay of the GI customers and t 2 is the mean delay 
of the Poissonian customers, then from (18) we have 

r 2 = -77^ — c + n a n (1 + MTi). 
m(1 - p) m(1 - P) 

Again, the second term on the right-hand side of this equation can be 
thought of as the added mean delay due to the presence of the GI 
customers. 

2.4 Measures of Variability of a Traffic Stream 

The delay suffered by an arbitrary customer in an input traffic-stream 
depends on how the customers' arrival epochs are distributed. Roughly, 
it can be said that the less "variation" in the arrival epochs, the better 
is the service received by the customers. We discuss two measures of 
this variation. 

Let A (£) be the interarrival time distribution of a traffic stream and 
Hi the jth moment of A(£): 

M , = r?dA®, j = 1,2, ••• . 

Jo 

We define V, the coefficient of variation of A(£), by 

V = 



Vl*a — Mi 



Mi 



This measure is dimensionless and depends only on the properties of 
the stream itself. 

Another measure which is used extensively in telephone traffic theory 
is denned with the aid of an infinite server system. The traffic stream 
is offered to an infinite number of servers with exponentially distributed 
service times. The ratio of the variance to the mean of the number of 
busy servers in statistical equilibrium is taken as a measure of the 
variation of the traffic stream. This number is called peakedness and 
is customarily denoted by Z. 
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We shall classify all traffic streams into two categories: smooth and 
peaked. A stream is smooth if V < 1 and peaked if 1 < V* The same 
dichotomy is effected by the inequalities Z < 1 and 1 < Z; however, 
Z is a function of the stream's intensity and thus is not as convenient 
in measuring variability properties. For example, if a(s) is the Laplace- 
Stieltjes transform of A(£), p is the stream's traffic intensity, and /z 1 
is the mean holding time, then it follows (Ref. 3, Chapter 3) that 

z - r=^ " '• (20) 

Using Jensen's inequality, 4 one can show that for a fixed X, the 
stream's peakedness (20) attains a minimum whenever A(£) is the 
one-point distribution 



A® = . 



0, (<\ 



In this case, 



at s = 



a(s) — e 



-(l/X). 



— n - x '" 



a(n) = e 

and with this substitution into (20) we see that Z can vary from \ 
(when p = oo ) all the way up to unity (when p = 0). In contrast, 
V = 0, independent of p. We conclude that Z is not a desirable measure 
of variation for smooth streams. We shall see later that for peaked 
traffic, Z turns out to be a good measure. Incidentally, the above argu- 
ment also shows that with exponential holding times the minimum 
possible value which peakedness can attain is £. 

2.5 Examples 

We now give two examples: a case with GI being smooth and another 
with GI being peaked. 

Example 1: GI Smooth. Let {X,} be a sequence of independent, expo- 
nentially distributed random variables with corresponding means { 6, } . 
Let the interarrival time distribution of the GI stream be given by the 

* A Poissonian stream (V = 1) is considered the norm to which the relative 
properties of smoothness and peakedness are compared. 



1316 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1972 

distribution of the sum 

Y = E x, , 

J=l 
where n > 1. Then the GI stream is smooth since 



In particular, if all the means are the same, then Y has a gamma distri- 
bution and GI is E n , an Erlangian stream. For this case, we have 

Vn 

which, moreover, is also the minimum value of (21). For our numerical 
example, we take n = 2. 
From (10) we have 

•« ■ tfeT 

and hence (9) becomes 

V d fephj) - (1 - i*)ff(g) 

with the derivative being evaluated at z = 2v. This follows from (14). 
Carrying out the differentiation, collecting terms, and solving for 
p(u), we obtain 

M = fg , (22) 

where 
R(u) = ia*{\itH l u a - [(X - n - 2>/)# 1 - tf 2 ]u + m^i}, 
T(li) = pV - (4<r 2 + p 2 + 2p + 4p«r)u 2 + (2p + 4<r + l)w - 1, 
X 



<;■ = — 



and 



H, = H'(2v), H 2 = H(2v). 
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It is easy to verify the following inequalities: 

T(0) < 0, T(l) > 0, T(\/ P ) < 0. 

Since p < 1 and T(u) is a polynomial of third degree, these inequalities 
imply that T(u) has three positive roots with one root being less than 
unity. Hence, if w, , u 2 , u 3 are the roots of T(u), then < u t < 1 < 
u 2 < u a . Since p(u) is a holomorphic function of u f or | u | ^ 1, the 
root Ui must also be a root of R(u). Using this requirement and the 
normalization condition p(l) = 1 we can find H t and H 2 . 
Omitting the intervening algebra we see that 

v(u ) = (1 ~ PU>v)K 

P{U) (1 - «,*)(! - o> 3 u) ' (23) 



where 



and 



£ _ (1 — (o 2 )(l — o> 3 ) 

(1 - pw x ) 



1 1 

«2 = , C0 3 = 

U 2 U 3 

Note that p , the probability of not being delayed at all, is given by K. 
Recall that for the Poissonian customers this quantity was given by 
q = 1 — p — <t. 
Expanding (23) in powers of u, we obtain 

P(u) = 7 v 2 [("2 - pwi)c4 - (« 3 - puM]u k , 

(w 2 — 0i s ) k=0 

and hence, the distribution [p f ] is given by 

Pi = -, r (w 2 — puM — ("a — pWi)w 3 , j = 0, 1, 2, • • • . 

The mean of this distribution can be computed by differentiation. 
We have 

0) 2 . (t>3 pUi 



V'iX) = 7-==- + 



1 — 0> 2 1 — C0 3 1 — pW t 

The corresponding quantities for the Poissonian stream of customers 
can easily be obtained using differentiation and relation (16). 

If D(t) is the delay distribution for the GI customers, then, as we 
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already have seen, its Laplace-Stieltjes transform is given by 

From (23) we have 

, M K(s + p)[s + (1 - pvM 

d{s) ~ [8 + (l -«0m][« + (i -« 8 )m]' 

Inverting this transform, we obtain the delay distribution for the GI 
customers: 



where 



D(t) = 1 - A 2 e- (1 - U ' , "« + A*" 11 — >", (24) 

j _ o> 2 (l — co 3 )(a>2 — pMi) 
(1 — pwi)(w 2 — co 3 ) 

. _ co 3 (l — fa> 2 )(<03 ~ A"l) 

(1 — pWi)(co 2 — co 3 ) 

Using (19) we can obtain E(t), the delay distribution for the Poissonian 
customers. Performing the indicated integration, we have 



E{t) = 1 - pe- {1 - p> " - a 1 + 



P — 0) 2 p — u 3 



~L-u-p>M< 



<rA 2 -(i-a.)m _ 3 g-' 1 -"''"' 1 (25) 

P — w 2 p — <o 3 

Figure 1 shows complementary delay distributions for the GI customers 
and the Poissonian customers in the E 2 + M/M/l model. Note that 
E 2 customers receive better service since they arrive in a smoother 
stream (V = l/vl2). For the Poissonian stream we, of course, have 
V = 1. While significant at low traffic intensities, this advantage 
diminishes as the traffic intensity increases. 

Example 2: GI Peaked. We can generate a peaked traffic stream using 
the interrupted Poisson process. 5 For a given a and Z, the interarrival 
time distribution of such a stream is given by Ref. 5. 

A(Q = h(l - e" r,{ ) + fc a (l - e~ r '*), 
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where 



and 




12 3 4 

T IN MULTIPLES OF MEAN HOLDING TIME 



Fig. 1 — Delay in the E 2 + M/M/l queue. 



Ti = *{t +W + 7 + V(r +co + T ) 2 - 4tco}, 
T % - \\r + « + 7 - V(t +« +t) 2 - 4tw{, 
t - r. 



fc, = 



r, — u 



1*2 — J- tC\ j 



T = aZ + 3Z(Z - 1), 



» - 1; - l h 

The Laplace-Stieltjes transform of A(£) is given by 
«(«)- 



"M^l I /C2^2 



s + r, s + r 2 

Carrying out the same steps as in the previous example, we see that 
the probability generating function of [p t ] has the same form, that is, 
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(1 - C0 2 )(1 - Q> 3 )(1 - pU t u) 



where 



P(W) " (1 - w 2 u)(l - w 8 «)(l - m) ' 



1 1 

OJ 2 = — . 0)3 — 



U 2 ' U a 

and Ui , u 2 , u 3 are now the roots of the following polynomial: 

p{p + fcin + fc 2 r 2 )w 3 - [ly, + fcif, + fc 2 r 2 + (ri + r 2 )p + 2p + p 2 ]w 2 

+ (r, + r 2 + 2 P + l)w - 1. 

Hence, the delay distributions for the GI customers and the Poissonian 
customers have the same forms as those in Example 1, but with different 
values of Wi , o> 2 , and w 3 ■ 

Figure 2 shows the same information as Fig. 1 of the previous example. 
Note that the peaked stream with Z = 3.0 (V = 2.3) receives poorer 
service than the Poissonian stream. Again, this effect diminishes at 
higher traffic intensities. 

Figure 3 shows the effect of peakedness and smoothness on the 
quality of service received by all customers. Holding the total traffic 
intensity constant (p + a = 0.8) and comparing the results with the 
case when all the customers arrive in a Poissonian stream (GI = M), 




10 20 30 40 

T IN MULTIPLES OF MEAN HOLDING TIME 



Fig. 2— Delay in the GI + M/M/l queue when GI is an overflow stream. 
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10 20 30 40 

T IN MULTIPLES OF MEAN HOLDING TIME 



Fig. 3 — The effect of peakedness on delay in the GI + M/M/l queue. 



we see that the service provided to all the customers either improves 
or deteriorates according to whether the GI stream is smooth or peaked. 

Figure 4 shows the effect of mixing the traffic. Note that as the 
proportion of peaked traffic in the total offered traffic stream increases, 
the service deteriorates for everyone. 

We now clarify a remark made earlier about the suitability of the 
peakedness factor Z for measuring variability properties of a stream. 
It was found numerically that the relation between Z and V in a peaked 
stream is monotone and nearly independent of the stream's intensity. 
In fact, by direct computation from A(£), we have the two relations 



1 + 7 2 



and 



Z = 



V = <1 + 



V 2 - 1 



3(1 + V 1 ) + 2a 



m - 1) 



1 - 1/(3Z + a) 



Figure 5 shows this relation for two streams, one with traffic intensity 
of 1 erlang and the other with 100 erlangs. Since the inverse relation 
is also monotone and nearly independent of the traffic intensity, the 
peakedness seems to be a suitable measure of the variational properties 
of a peaked stream. 
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Fig. 4— Effect of mixture of input on delay in the GI + M/M/l queue. 
III. THE GI + M/M/ oo SYSTEM 

Let Y(t) be the number of customers in the system at time t. Since 
there are an infinite number of servers, Y(t) is equal to the number of 
busy servers at time t. Again, we see that { Y(t), t ^ 0} is a piecewise 
Markov process. The regeneration points are the arrival epochs of 
the GI customers; the distribution of the length of Markovian segments 
is given by A(£), the interarrival time distribution of the GI customers; 
and the regeneration matrix is given by (1). This is the same identifica- 
tion as that made in the GI + M/M/l queue. The difference here is 
the Markov process operating within the segment. 

The Markovian development of the process within the segments is 
governed by a birth-death process with birth rate X, the arrival in- 
tensity of the Poissonian stream, and death rate j/t when the process 
is in state j, where /T 1 is the mean service time. The transition functions 
of this process are given by the transient solution to the M/M/oo 
system (Ref . 2, page 24) : 



rM-ZdlJk'-™ 



[1 - rf0]"V"' + "<0, 



(26) 



where 



g(t) - 1 - e 



P = 
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Let {p,} be the stationary distribution of the Markov chain imbedded 
at points immediately preceding a GI arrival. This distribution satisfies 
the Chapman-Kolmogorov eqs. (3). If p(u) is the probability generating 
function of {p,}, then from (4) we have 



P(u) = f 't,PiG i+l (u,QdA(Q, 

Jo t=0 



where 



J=0 

But from Ref . 2 we have 

G<(u, Q = (1 - e" M£ + we-^V'"-"'" 1 *', 
and we see that p(w) satisfies the following integral equation: 

p(u) = e- p(1 - u) f r(u, Mr(u, ^-^^^ d A{Q, 
Jo 



where 



r(u, Q = 1- e-" k + ue-'t. 



(27) 



Now let Xi be the number of Poissonian customers and X 2 the 
number of GI customers in the system seen by an arbitrary GI arrival. 
Since the number of servers is infinite, X Y and X 2 are independent. 
If b(u) is the probability generating function of the distribution of X 2 , 
then 



p(u) = e 



-P(l-u) 



b(u). 



(28) 
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Fig. 5 — Relation between V and Z in an overflow stream for two different means. 
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Since the arriving GI customers are sampling a Markov process ac- 
cording to a renewal process, they will see the stationary distribution 
of the process 1 and, hence, X, has the Poisson distribution 



k 



PiXt-k) =e- p £y, (29) 

the stationary distribution of the number of customers in an M/M/ oo 
system. 

Substituting (28) into (27) we obtain an equation for b(u): 

b{u) = f r(u, ®b[r(u, £)] dA(£). (30) 

Jo 

In solving for b(u), it is convenient to consider the expansion of b(u) 
about u = 1, rather than about u = 0. Hence, we set 

b(u) = Ed,(w- 1)', 

»=o 

where 
If 

1=0 

is the expansion of b(u) about the origin, then the two sets of coefficients 
are mutually related through the equations 

h = t ft)(-i) fc "' 4 • ( 31 ) 

4 = 1 W' 

Differentiating (30) j times and setting u = 1, we get 
dj = djoj + df-ia,- » J — li 2, • ■ ■ , 
d = 1, 

where 



a, =a(j M ) = Ce-^dAiO. 
Jo 
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Solving for d, , we obtain 

di= flr=T' 3-1,2,.-- . 
Thus, we have 

p(u) = e -" (1 - u) £ chiu - 1)' 

1=0 
1=0 

and the distribution {p,-} is given by the convolution of {&,-} with (29): 

*=o 'C! 

Now 

g, = limP{r(0 = il, J = 0,1,2, ••• , 

the distribution of the number of busy servers seen by an arbitrary 
Poissonian arrival, can be found using the rate conservation principle. 
If v~ ' is the mean of A (£) and a = v/n, then equating the asymptotic 
rate of transition out of the set of states {0, 1, 2, • • • , j) to the rate 
into that set we get 

*Vi + PQi = (j + l)ft + i , j = 0, 1, 2, • • • . (32) 

If q(u) is the probability generating function of {#,}, then it follows 
from (32) that q(u) satisfies the differential equation 

q'(u) = pq(u) + <rp(u). 

This equation has the general solution 

q(u) = e'V / e-"p(Q d£. 

Substituting for p(£), carrying out the integration and using the nor- 
malization condition q(l) = 1 to determine the integration constant, 
we obtain 

Since q(u) is the product of two probability generating functions, the 
distribution j g, } can now easily be determined by convolution. Another 
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way of obtaining {g,j is by solving (32), i.e., 

ft-^ + tE^i i= 1,2, -•■ , 



yTo ]\ hd P 
The means of the two distributions are given by 






2 hi = 9'(1) = P + a 
j-i 

We make two observations. First, the mean of fg,} is independent of 
the form of the interarrival time distribution of the GI customers; 
it depends only on v~ l , the mean interarrival time. Second, using 
Jensen's inequality one can show that the mean of (p,} is minimized 
whenever the interarrival times of the GI customers are constant. 
For this special case we have 



and hence 



since 



for a > 0. 



e ua - 1 
ff'(l) - P + a, 

P'(D < <z'U) 



1 + - < e 1/0 
a 
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